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Abstract 

The  development  of  fuel  cell  cars  can  play  an  important  role  in  resolving  transport  problems,  due  to  the  high  environmental  compatibility 
and  high  efficiency  of  this  kind  of  vehicle.  Among  the  different  types  of  fuel  cells,  proton-exchange  membrane  fuel  cells  (PEMFCs)  are 
considered  the  best  solution  for  automotive  applications  at  the  moment. 

In  this  work,  constructive  criteria  are  discussed  with  the  aim  of  obtaining  a  power  generation  module  adaptable  to  a  wide  range  of  cars. 

A  particular  problem  in  accomplishing  the  overall  project  is  represented  by  the  definition  of  the  compressor  system  for  air  feeding.  In  this 
work,  the  design  approach  to  the  problem  will  be  delineated:  some  options  are  reviewed  and  the  best  solution  is  analysed. 

The  transient  response  of  the  system  (fuel  cell  and  compressor)  is  investigated  in  order  to  optimise  the  start-up  running  through  a  model  of  a 
fuel  cell  stack  and  a  compressor  simulation. 

The  model  and  its  results  are  proposed  as  a  work  procedure  to  solve  the  problem,  by  varying  external  conditions:  in  fact,  to  perform  the 
system  start-up  under  stable  conditions,  the  air  relative  humidity  and  temperature  must  be  maintained  in  a  proper  range  of  values. 

The  approach  here  presented  has  been  utilised  for  the  definition  of  the  characteristics  of  the  power  module  and  layout  of  a  middle-size 
hybrid  city  bus  in  the  framework  of  a  project  promoted  by  the  European  Union.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

The  increase  in  pollutant  emission  levels  and  the  parallel 
development  of  environmental  awareness  have  forced  the 
automotive  industry  to  study  innovative  solutions  to  improve 
the  quality  of  the  urban  environment.  The  technical  problem 
was  approached,  above  all,  by  US  researchers  in  the  early 
1980s  and  the  solutions  realised  consisted  of  the  control  of 
the  emissions  of  traditional  vehicles  (catalysers  and  parti¬ 
culate  recovery  systems). 

This  approach  achieved  a  notable  decrease  in  toxins, 
pollutants  and  irritant  components,  such  as  CO,  NOx  and 
hydrocarbons.  Following  US  laws,  the  European  directives, 
which  in  fact  only  impose  the  catalytic  conversion  of  emis¬ 
sions,  also  support  the  studies  of  non-pollutant  vehicles. 

At  the  moment,  electric  vehicles  have  the  largest  share  of 
the  low-pollutant  vehicle  market  but  their  autonomy  is  still 
not  sufficient  for  normal  use  and  the  area  required  for  the 
batteries  considerably  reduces  the  usable  space. 

With  the  aim  of  obtaining  vehicles  with  the  intrinsic 
possibility  of  zero  emissions,  new  solutions  of  the  problem 
were  examined  and  tested  in  the  early  1990s. 
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One  of  the  most  promising  technologies  for  obtaining  the 
desired  result  is  represented  by  fuel  cells,  where  electrical 
energy  is  obtained  by  direct  conversion  of  combustible. 
Industrial  solutions  using  this  electrochemical  device  in 
the  automotive  industry  are  the  proton-exchange  membrane 
fuel  cells  (PEMFCs),  producing  electrical  energy  at  low 
working  temperature  with  high  efficiency  [1]. 

Furthermore,  these  systems  offer  the  best  solution  for 
reducing  pollution  to  zero  in  city  centres:  the  fuel,  usually 
hydrogen,  is  burned  producing  only  water.  This  solution  by¬ 
pass  NOx  production  problems,  typically  involved  in  tradi¬ 
tional  engines  powered  by  hydrogen. 

The  definition  of  the  general  layout  criteria  of  a  fuel  cell 
vehicle  presented  in  this  work  follows  a  method  that  can 
offer  great  advantages  in  terms  of  assembling  simplicity 
(possibility  of  endowing  an  existing  car  with  a  power 
module)  and  modularity  but  it  shows  some  start-up  pro¬ 
blems.  In  fact,  following  this  approach,  the  power  generation 
is  influenced  by  the  performances  of  the  motor  of  the 
compressor.  In  addition,  the  pressure  inside  the  stack  is 
furnished  by  the  compressor  and  its  motor  is  fed  by  the 
energy  produced  by  the  stack. 

The  need  to  examine  the  transient  response  of  the  fuel  cell 
stack  coupled  to  the  air  compressor  has  been  resolved  with  a 
model  of  the  transient  behaviour  of  the  system,  developed  on 
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Nomenclature 

a  relative  humidity,  defined  by  Eq.  (A. 3) 

Cp  mass  specific  heat  (J  kg  1  K1) 

D  diffusivity  (m2  s  1 ) 

'y  _  | 

D2  diffusivity  as  a  function  of  /l  (m  s-  ) 

F  Faraday  constant  (C  mol  :) 

AG  Gibbs  free  energy  change  of  the  overall 
reaction  (J  mol  x) 

AGst  standard  AG  of  the  overall  reaction  (J  mol  x) 
h  heat  transfer  coefficient  (W  m  2  K-1) 

AH  enthalpy  change  of  the  overall  reaction 
(J  mol-1) 

I  electrical  current  density  (A  m  2) 

70  exchange  current  density  (A  m  2) 

<y 

k  hydraulic  permeability  of  the  membrane  (m  ) 

kh  Boltzmann  constant  (J  K-1  per  molecule) 

m  membrane  thickness  (m) 

M  equivalent  weight  of  membrane,  i.e.  weight  of 

membrane  per  mole  of  S03-  groups  (g  mol-1) 
Mt  total  weight  of  stack  (kg) 

n  number  of  moles  (mol) 

ne  number  of  electrons  transferred  in  the  reaction 

Pi  partial  pressure  (Pa) 

P  total  pressure  (Pa) 

R  gas  constant  (J  mol-1  K-1) 

sair  contact  area  air  per  cell  (m2) 

sest  external  surface  (m2) 

S  total  internal  cell  surface  (m2) 

t  time  (s) 

T  temperature  (K) 

V  potential  (V) 

Vs  inner  stack  volume  occupied  by  air  (m3) 

v  mole  fraction 

z  co-ordinate  (m) 

Greek  letters 

oc  parameter  defined  by  Eq.  (A. 5)  (mol  m  2  s-1) 

P  charge  transfer  coefficient 

p  parameter  defined  by  Eq.  (A. 6)  (mol  m  s  ) 

X  ratio  net  water  flux  through  the  membrane/ 

water  produced  by  the  electrochemical  reaction 

8  parameter  defined  by  Eq.  (A.  8)  (mol  m  2  s-1) 

<P  coefficient  in  Eq.  (10) 

y  parameter  defined  by  Eq.  (A. 7)  (mol  m  1  s-1) 

yi,  y2,  coefficients  in  Eq.  (13) 

73,  74 

r)  overpotential  (V) 

(pi,  cp2,  coefficients  in  Eq.  (A.4) 

(Pl,  (P4, 

(P5,  (P6 

X  ratio  H20/S03-  in  membrane 

X0  value  of  X  at  the  cathodic  side 

X\  value  of  X  at  the  anodic  side 

9  diameter  of  water  molecule  (m) 


p  membrane  density  (kg  m  3) 

o  resistivity  (Q  1  m-1) 

t  vapour  pressure  of  water  (Pa) 

d  stoichiometric  coefficient 

Q  electrical  resistance,  Q  m 

£  parameter  defined  by  Eq.  (A. 9)  (m-1) 

i fa,  \jj2,  coefficients  in  Eq.  (A. 2) 

<A3,  *>4 

C  electro-osmotic  drag  coefficient 

Subscripts 
a  anode 

c  cathode 

air  air 

est  external  stack  surface 

exp  experimental 

ext  external  to  the  stack 

H20  water 

02  oxygen 

stack  stack 

u  total  number  of  moles 

343  evaluated  at  temperature  of  343  K 

Superscripts 

0  initial  conditions 

c  cell 

el  electrode 

Is  losses 

memb  membrane 
t  thermodynamic 

the  basis  of  experience  acquired  through  the  participation 
of  the  group  of  University  of  Genoa  in  the  European  project 
for  the  development  of  a  middle-size  city  bus. 

In  general,  at  high  values  of  temperature  and  relative 
humidity  the  power  furnished  by  the  fuel  cell  stack  (strongly 
dependent  on  these  variables)  is  sufficient  to  move  the  motor 
of  the  compressor;  moreover,  the  transient  model  of  the 
entire  system  allowed  us  to  analyse  some  running  options 
and  to  determine  the  best  solution  for  the  start-up  of  the 
“power  module  system”. 

2.  General  description  of  a  fuel  cell  vehicle 

The  solution  here  proposed  for  the  realisation  of  a  fuel  cell 
vehicle  can  be  considered  an  innovative  approach  to  the 
problem. 

In  fact,  hybrid  fuel  cell  vehicles  [2]  are  usually  realised  by 
coupling  an  electric  power  generator  to  a  battery  storage 
system.  The  surplus  energy  is  stored  and  utilised  when  the 
electric  motor  needs  peak  power. 

In  the  solution  proposed  here,  the  power  generator  is 
utilised  only  to  charge  the  batteries  which  feed  the  engine 
of  the  bus.  In  addition,  two  alternative  solutions  for  the 
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energy  supply  to  the  motor  of  the  compressor  are  reviewed: 
Fig.  1A  represents  the  “traditional”  approach,  feeding  the 
motor  from  the  batteries  [2];  Fig.  IB  represents  the  novel 
approach  where  the  motor  of  the  compressor  is  fed  by  the 
power  generator. 

A  hybrid  fuel  cell  vehicle  designed  following  the  latter 
approach  can  be  divided  into  two  distinct  sections:  the  first 
section  consists  of  the  chassis  with  steering,  batteries, 
electric  motor,  braking  system  and  electronic  control  centre 
for  all  equipment  and  the  second  one  is  represented  by  the 
energy  generation  system,  consisting  of  one  or  more  fuel  cell 
stacks  and  auxiliaries  (air  compressor,  hydrogen  or  other 
fuel  feeder,  heat  exchangers,  dc/dc  converter  and  a  max¬ 
imum  power  point  tracker  to  feed  the  batteries).  The  scheme 
of  the  bus  reported  in  Fig.  2  highlights  the  innovative 
approach. 

The  modularity  of  the  system  permits  the  construction  of  a 
pack  independent  of  battery  characteristics,  in  agreement 
with  the  typical  modular  concept  of  fuel  cells.  Following  this 
criterion,  it  is  possible  to  use  any  chassis  as  a  base,  making 
sure  to  put  the  electrical  generation  system  in  an  autonomous 
chassis  connectable  to  the  principal  one. 


2.1.  Fuel  cell  power  generator 

PEMFCs  are  actually  considered  the  best  solution  for 
automotive  applications  [1].  This  is  due  to  criteria  related  to 
the  flexibility  of  their  use  (starting  is  very  simple,  the  device 
is  of  small  dimensions,  efficiency  is  very  good)  and  the  low 
working  temperature,  which  is  important  for  the  intrinsic 
safety  of  the  vehicle. 

These  cells  are  formed  of  a  perfluo sulfonic  membrane  as 
electrolyte  coupled  with  two  electrodes  (anode  and  cathode) 
built  utilising  a  substrate  of  carbon-PTFE  coated  with 
platinum  as  electrocatalyst. 

Air  is  fed  to  the  cathodic  compartment  while  hydrogen  is 
fed  to  the  anodic  one  and  the  electrolyte  performs  both  the 
functions  of  transferring  H+  from  the  anode  to  the  cathode 
and  reactant  separation. 

Electrochemical  reactions  involved  are: 

H2  — >  2H+  +  2e“  (anode) 

2H+  +  ^02  +  2e“  — >•  H20  (cathode) 

H2  +  ^02  — >•  H20  (overall  reaction) 
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Stacks  are  obtained  by  connecting  a  number  of  cells  in 
series  while  gas  distributors  (the  bipolar  plates)  accomplish 
the  function  of  feeding  the  fuel  to  the  anodic  side  and  air 
(or  oxygen)  to  the  cathodic  one. 

Using  common  relationships  between  the  power  continu¬ 
ously  supplied  by  the  generation  module  and  that  adsorbed 
by  the  vehicle  and  paying  particular  attention  to  one  of  the 
possible  urban  cycles — in  this  case  represented  by  Cole’s 
relationship  [3] — the  peak  power  requested  by  a  middle-size 
city  bus  can  be  considered  about  100  kW,  while  the  power 
consumed  by  the  auxiliary  systems,  including  all  the  losses, 
can  be  estimated  at  about  10  kW.  So,  the  battery  design  is 
based  on  the  maximum  power  needed  by  the  whole  system 
(bus  and  auxiliaries).  As  a  consequence,  to  guarantee  con¬ 
tinuous  bus  running,  the  installed  power  furnished  by  the 
fuel  cells  must  be  about  50  kW  [4]. 

In  this  case,  the  choice  of  the  low  power  installed  in  the  bus 
is  due  to  economical  reasons  in  the  development  of  the 
prototype.  In  addition,  due  to  the  characteristics  of  the  mod¬ 
ularity  of  this  solution,  the  results  obtained  will  be  adaptable 
to  a  wide  number  of  automotive  realisations  (cars  and  buses). 

2.2.  Reactant  feeding  system 

The  supply  of  fuel  (H2)  to  the  cells  is  usually  a  simple 
operation  from  the  point  of  view  of  the  regulation  of  the 
pressure  and  flow  rate.  In  fact,  the  gas  is  contained  in  one  or 
more  cylinders  on  the  roof  of  the  vehicle  and  through  valve 
control  it  is  possible  to  obtain  the  characteristic  desired. 

On  the  contrary,  the  supply  of  air  is  not  so  straightforward: 
it  is  necessary  to  feed  a  constant  air  flow  rate  with  a  very  high 
purity  level  to  avoid  poisoning  the  fuel  cells. 

So,  the  possible  options  connected  with  the  choice  of  the 
air  compressor  system  refer  to  working  pressure,  air  purity 
and  system  layout. 


Literature  data  on  the  fuel  cell  power  supply  [5]  show  the 
best  performances  for  high  air  pressure,  but  to  get  a  simple 
constructive  solution  and  management,  low  working  pressure 
has  been  chosen:  the  decrease  in  electrical  performances  is 
compensated  by  the  simplicity  of  the  system  and  the  contain¬ 
ment  of  costs. 

With  regard  to  the  problem  of  air  purity,  the  solution 
that  guarantees  the  complete  absence  of  dragged  oil  or 
metallic  shards  is  represented  only  by  an  oil-free  screw 
compressor.  This  solution  offers  all  the  advantages  of  the 
rotative  compressors  without  wear,  drag  and  eventual  abra¬ 
sion  of  moving  parts.  The  lack  of  wear  is  achieved  thanks  to 
two  rotors  with  low  but  sufficient  tolerances  to  guarantee 
no  contact.  High  rotational  velocity  (about  30,000  rpm) 
provides  the  chosen  working  pressure  in  only  one  stage. 
The  characteristic  curves  of  the  compressor  installed  in  the 
bus  are  reported  in  Fig.  3:  the  solid  line  represents  the  power 
needed  by  the  compressor  as  a  function  of  its  rotation  speed; 
the  broken  line  the  air  flow  rate  and  the  dotted  line  the 
adiabatic  efficiency.  The  maximum  adiabatic  efficiency  is 
reached  when  the  compressor  works  at  60%  of  the  nominal 
power. 

In  the  layout  of  the  air  feeding  system  described  earlier 
(Fig.  IB),  all  the  devices  present  in  the  power  module  can 
only  take  energy  from  the  electrochemical  generator. 

3.  Simulation  of  the  power  module 

A  wide  number  of  works  [6-12]  demonstrate  that  the 
starting  of  PEMFCs  is  instantaneous,  but  practically  no 
information  is  available  on  the  start-up  problems  of  the 
entire  power  module  system  (stack  +  compressor  +  motor). 
To  analyse  the  transient  behaviour,  study  the  influence  of 
environmental  conditions  and  determine  the  best  start-up  path 
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without  operating  fuel  cell  devices  under  limit  conditions,  a 
simulation  model  of  the  power  module  system  was  set  up. 

The  simulation  of  a  power  module  includes  both  the 
model  of  the  fuel  cell  stack  and  the  experimental  trend  of 
the  compressor  and  of  the  motor  system.  This  “mixed” 
approach  is  motivated  by  the  need  to  predict  the  stack 
operative  conditions  that  will  avoid  damaging  drying  or 
flooding  of  the  membrane. 

3.1.  Stack  simulation 

Theoretical  approaches  [7,8,1 1-14]  have  often  been  devel¬ 
oped  under  steady  state  conditions  and  in  general  are  not  able 
to  describe  fuel  cells  coupled  with  an  external  utilisation. 

This  work  presents  a  dynamic  study  of  an  electrochemical 
device  which  includes  a  macroscopic  simplified  fuel  cell 
model  coupled  with  a  detailed  description  of  the  membrane. 
This  tool  furnishes,  in  reasonable  time,  a  satisfactory 
description  of  the  stack  behaviour  comparable  with  experi¬ 
mental  literature  data  and  more  complex  approaches. 

The  main  assumptions  of  the  model  are: 


The  time -dependent  equations  of  the  model  are  written  in 
a  macroscopic  form  and  in  order  to  have  a  fast  numerical 
integration  including  mass  balances  inside  the  cell  and  the 
energy  balance  of  the  solid. 

Cell  voltage  is  calculated  taking  into  account  activation 
overpotential,  Ohmic  effects  and  diffusive  limitations.  Then, 
the  results  of  a  single  fuel  cell  have  been  extended  to  an 
entire  stack  considering  this  formed  by  a  number  of  cells 
with  the  same  properties. 

The  time-dependent  that  assume  the  feeding  pressure  of 
the  hydrogen  constant  for  every  operative  condition,  are 
reported  in  the  following: 


•  for  water 


<9«h2o 

dt 


- S  if  a<  1 

neF 


(1) 


or 


dnn2o  _ 
dt 


if  a  >  1 


(2) 


•  uniformity  of  the  stack  temperature; 

•  hydrogen  exiting  from  its  bottle  instantaneously  reaches 
the  stack  temperature  (this  is  due  to  the  low  hydrogen 
specific  heat); 

•  momentum  balance  is  substituted  by  the  experimental 
data  obtained  from  a  screw  compressor  coupled  with  a 
fuel  cell  stack  simulator; 

•  operating  conditions  far  from  diffusion  limiting  ones. 


with  initial  conditions: 


for  t  =  0,  n\0 

•  for  oxygen 
dno2  IS 


/4o  0 

p0  u 


dt 


2nPF 


(3) 

(4) 
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with  initial  conditions: 


P°V 

for  t  =  0,  nQ  = - -  and  nl\  = 

■  U  RT  Ui  U  U2 


(5) 


Some  considerations  are  necessary  for  the  interpretation 
of  the  energy  balance  of  the  solid.  The  evaluation  of  start-up 
effective  strategies  are  evaluated  under  safe  conditions,  i.e. 
under  the  worst  conditions  for  the  whole  system:  absence  of 
external  coolant/heating,  heating  by  feeder  gas  only  during 
the  steady  state  compression  phases,  absence  of  forced 
humidification. 

For  these  reasons,  the  energy  balance  is  interpreted  as  a 
time-dependent  equation  that  relates  stack  temperature  var¬ 
iations  to  the  heat  losses  and  the  exchanges  with  the  air  flow 
and  the  environment: 


Mt(Cp)  stack 


^T^stack 


dt 


^aiAair  (Tiir  —  Tstack) 


+  ^est^est(7ext  — T’stack)-  VC^jlS 

(6) 


The  water  condensation  in  the  energy  balance  has  been 
neglected  because  the  term  of  condensation  enthalpy  is 
negligible  if  compared  with  the  others.  Heat  transfer  coeffi¬ 
cients  are  evaluated  as  that  of  the  forced  air  for  the  inner  side 
and  the  natural  convection  for  the  external  one,  following 
Amphlett  et  al.  [6].  For  the  same  reasons  as  discussed  above 
no  coolant  or  heater  are  present. 

The  calculation  of  the  Nernst  potential  (V): 


following  Bemardi  and  Werbrugge  [13]  and  Springer  et  al. 
[14]: 


RT  in  (J  P  \ 

w  Uw 


(9) 


where  the  ratio  P/po2  includes  the  concentration  overpo¬ 
tential  and  7//0  the  activation  one. 

Exchange  current  density  (I0)  is  evaluated  following 
Anderson  et  al.  [15]: 


T-343 

1 T  —  ^0 1 343  exp 

(p - 

343T 

(10) 


Then,  the  evaluation  of  the  overall  Ohmic  effects  in  the 
membrane  using  a  monodimensional  model  follows 
Springer  et  al.  [14,16]:  Ohmic  losses  are  expressed  as  a 
function  of  the  internal  membrane  humidity  {X),  this  related 
to  the  total  current  and  the  water  contained  in  the  active  sites 
of  the  membrane  (S03  for  a  perfluosulfonic  membrane). 

The  membrane  Ohmic  losses  are: 


^memb  =  Qj 


(ID 


where  the  Ohmic  membrane  resistance  is  calculated  along 
the  membrane,  by  the  integration  of  the  conductivity: 


(12) 


Eq.  (12)  is  utilised  by  considering  the  conductivity  taken 
from  experimental  evidence  (Appendix  A),  as  a  function  of 
the  ratio  between  the  water  and  the  S03  in  the  membrane 
and  the  mean  temperature: 


a  =  (yxX  -  y2)  exp 


f3 


(13) 


is  the  basis  for  the  evaluation  of  the  energetic  performances 
of  the  single  cell. 

In  fact,  in  the  PEMFCs  the  open  circuit  voltage  does 
not  correspond  to  the  Nernst  voltage:  this  is  due  to  parasi¬ 
tic  phenomena,  such  as  undesirable  electrochemical  reac¬ 
tions  at  the  cathode  which  produce  an  overpotential  that 
is  evaluated  as  a  constant  value  in  accord  with  Rho  et  al. 
[11]: 

rjh  =  0.2V  (8) 


where  X  is  conditioned  by  the  electro-osmotic  drag,  current 
density  and  pressure  gradient  between  the  two  electrodic 
sides: 

#  =  x(i—2L _ P*”'-1'  )  _  i——  (14) 

dz  V  neF  pDx  khT  m  )  A  neF  pDx 

Taking  into  consideration  that  diffusivity  (Dp  can  be  con¬ 
sidered  constant  when  no  limiting  current  is  reached  [14], 
Eq.  (14)  can  be  integrated  analytically  to  obtain  an  explicit 
form  of  Ohmic  resistance  (Appendix  A): 


exp  [y(TW(£,m)]{£m-]n  [(y^i  ~f2)(ex P  (<»  ~  1)1  +in  V(^m)gW}} 

C[exp(<»g(2o)  -  X\yx  +  y2] 


The  effective  voltage  can  be  calculated  considering  the  main 
losses. 

Losses  occurring  at  the  cathode  are  mainly  due  to  activa¬ 
tion  effect  and  they  are  of  one  order  of  magnitude  higher 
than  the  overall  anodic  overpotential  [12]:  so,  anodic  over¬ 
potential  will  be  neglected.  Cathodic  activation  overpoten¬ 
tial  and  concentration  overpotentials  have  been  calculated 


As  evidenced  in  Eq.  (15)  Ohmic  resistance  is  conditioned  by 
the  relative  humidity  of  the  reactants  (20  and  Xp. 

Then,  the  cell  voltage  is  expressed  by: 

Vc  =  Vl-  4s  -  if  -  timemb  =  /(/ )  (16) 

The  model  is  integrated  following  two  paths.  In  the 
former,  able  only  to  give  a  static  interpretation  of  the  stack 
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Table  1 

Parameters  of  the  simulation  model 


Parameter 

Value 

AA43 

0.9  A  m~2 

P° 

1  x  105  Pa 

(p 

35000 

7i 

0.58  Q^1  irT1 

72 

0.50  Q-1  irT1 

73 

1268  K 

74 

303  K 

<Pi 

5.1  x  10-10m2s-1 

(p2 

-6.6  x  10~n  m2  s_1 

(P3 

5.3  x  10~12  m2  s_1 

<P4 

-1.3  x  10~13  m2  s-1 

(P5 

2416  K 

<P6 

303  K 

<Al 

0.043 

•/c 

17.81 

' h 

-40 

• h 

36 

?/s 

0.2  V 

p 

2.9  x  10  1  kPas"1 

power  as  a  function  of  external  conditions  (temperature  and 
relative  humidity),  the  method  of  finite  differences  has  been 
adopted.  In  the  latter,  where  it  was  necessary  to  interpret  the 
transient  response,  finite  differences  were  coupled  with  a 
relaxation  method  for  time-dependent  equations. 

The  code  is  written  in  FORTRAN  language  and  the 
resulting  software  runs  on  a  personal  computer  in  less  than 
a  minute. 

The  characteristic  curve  of  a  PEM  commercial  stack 
which  consists  of  100  cells  (surface  about  0.25  m2  per  cell) 
is  evaluated  for  average  cell  temperature  (T  =  343  K,  pair  = 
3.5  x  105  Pa,  pn2  =  4  x  105  Pa  and  a  =  0.99)  utilising  the 


parameters  reported  in  Table  1 ;  then  the  result  referred  to  the 
single  cell  is  reported  in  Fig.  4  to  be  compared  with 
experimental  data  reported  in  literature  [5,12]  obtained 
under  the  same  conditions. 

The  model  contains  simplifications  that  produce  a  5% 
difference  in  the  electric  response  of  the  stack  compared  to 
the  results  of  other  more  elaborate  models  but  the  differ¬ 
ences  obtained  by  comparing  the  proposed  model  with  some 
experimental  results  [8,11,12]  are  of  the  same  magnitude 
obtained  by  applying  the  other  models  to  these  data 
[6-8,11,12,14,16].  This  result  confirm  the  validity  of  the 
model  presented  in  this  work. 

The  simplifications  provide  an  analytic  solution  for  the 
Ohmic  overpotential  that  significantly  reduces  the  computa¬ 
tion  time  which  decreases  dramatically  from  10  min  for  the 
calculation  of  only  a  point  of  electrical  response  ( I-V)  to  a 
few  seconds  for  the  evaluation  of  the  entire  electrical  curve. 

In  addition,  the  electric  response  of  the  stack  has  been 
calculated  in  different  conditions  of  external  temperature 
and  air  relative  humidity:  as  can  be  seen  in  Fig.  5,  the  power 
generation  is  strongly  dependent  by  these  two  variables. 

Power  values  are  represented  as  a  percentage,  considering 
100%  that  is  furnished  at  T  =  343  K  and  100%  of  relative 
humidity  (about  30  kW). 

The  velocity  of  execution  of  the  model  permit  its  use  to 
follow  the  transient  behaviour  of  a  complex  system  coupling 
the  obtained  results  with  the  experimental  trends  of  the 
compressor. 

3.2.  Compressor  and  motor  experimental  trends 

To  accomplish  the  simulation  of  the  power  module  and  to 
analyse  the  start-up  of  the  whole  system,  the  power  needed 


current  density  [A  m'2] 


Fig.  4.  Comparison  between  results  and  experimental  literature  data  obtained  under  the  same  operating  conditions. 
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Fig.  5.  Stack  model  output:  results  of  power  generation  as  a  function  of 
temperature  and  relative  humidity. 


were  obtained  in  graphic  format  as  a  function  of  the  rotations 
of  the  motor. 

All  device  characteristics  that  were  measured  when  the 
stack  valves  were  opened,  are  reported  in  Fig.  6. 

So,  the  power  requested  by  the  motor  is  composed  with 
that  produced  by  the  fuel  cell  stack  to  individuate  the 
possible  start-up  strategies. 

The  values  of  air  temperature  at  the  outlet  of  the  com¬ 
pressor  are  utilised  for  the  evaluation  of  the  stack  heating 
during  transients,  so  the  pressure  values  in  the  stack  are 
available  and  when  the  stack  is  connected  to  the  compressor, 
these  data  permit  the  calculation  of  the  pressure  in  the  stack 
introducing  the  experimental  pressure  data  (Pexp): 


P{t)  —  P  exp  (t)  + 


fr  fdn02  <9nH2o\ 

Jo  V  dt  +  dt  ) 


(17) 


by  the  motor  of  the  compressor  must  be  known.  These  data 
have  been  collected  through  participation  in  the  European 
organisation  for  the  construction  of  a  middle-size  city  bus. 

Experimental  data  were  obtained  in  an  industrial  labora¬ 
tory  using  a  procedure  that  involved  the  connection  of  a 
screw  compressor,  fed  by  a  motor,  to  a  simulator  of  exhaust 
loss  (a  cylinder,  with  a  calculated  number  of  tubes  passing 
through  it,  through  which  air  was  forced)  designed  to  the 
specifications  of  a  30  kW  fuel  cell  stack. 

The  motor  was  fed  by  an  external  connection.  Some 
sensors  were  positioned  to  measure:  power  adsorbed  by 
the  motor,  pressure  in  the  mean  and  at  the  end  of  the  stack, 
temperature  and  air  flow  rate  at  the  compressor  escape 
flange. 

Measures  were  recorded  by  connecting  all  the  sensors  to 
an  oscilloscope  Tektronix  model  TDS  420A  and  all  data 


In  this  way  the  transient  stack  model  is  coupled  with  the 
motor  +  compressor  response  to  obtain  the  transient  beha¬ 
viour  of  the  whole  power  generator  module. 

4.  Discussion  of  start-up  criteria 

Four  start-up  strategies  have  been  investigated  utilising 
the  model  described  above. 

In  the  first  it  is  considered  that  during  the  motor  start-up 
(from  0  to  235  ms)  air  contained  in  the  stack  is  not  fed  by  the 
compressor  flow.  Air  pressure  and  oxygen  concentration 
vary  as  a  function  of  the  total  energy  furnished  to  the  motor 
(including  losses).  Steady  state  is  reached  when  the  time  is 
greater  than  235  ms:  only  after  this  time,  do  pressure  and 
oxygen  concentration  become  equal  to  the  equilibrium 


Fig.  6.  Motor-compressor  start-up  behaviour. 
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Fig.  7.  Temperature  and  relative  humidity  start-up  limits  for  strategies  no.  1  and  2.  Above  the  two  lines,  start-up  can  be  realised.  Below  the  two  lines,  start-up 
cannot  be  reached  without  humidification  or  without  heating  the  stack. 


conditions,  considering  the  experimental  values  for  the  stack  In  all  the  three  former  solutions,  it  is  supposed  that  the 

out  of  line  with  the  internal  consumption.  humidifier  is  absent:  results  are  conservative  and  discrimi- 

The  second  strategy  considers  that  air  pressure  on  the  nate  external  conditions, 
stack  follows  the  experimental  values;  oxygen  concentration  The  fourth  solution  is  analogous  to  the  third  but  relative 

is  evaluated  as  above.  humidity  is  maintained  >50%,  utilising  an  external  humi- 

The  start-up  can  be  realised  in  safe  conditions  only  when  difier.  The  results  of  the  two  latter  strategies  are  reported  in 

the  power  furnished  by  the  cell  is  greater  that  those  utilised  Fig.  9,  realised  in  the  same  way  as  Fig.  7. 
by  the  motor  of  the  compressor  in  all  the  early  stages.  Fig.  7  The  simulation  shows  that  the  novel  approach  of  Fig.  IB 

highlights  the  start-up  limits  of  these  two  solutions.  can  be  adopted  and  through  the  simulation  of  a  stack  in  the 

In  the  third  simulation  strategy,  a  multi-step  heating  has  absence  of  a  humidifier,  the  operative  limits  of  the  solutions 

been  considered:  the  motor  reaches  the  maximum  power  after  are  identified  (i.e.  external  temperature  and  air  relative 

the  path  described  in  Fig .  8 ;  during  the  steady  state ,  the  stack  is  humidity) . 

heated  by  reactions  and  gas  heat.  In  the  transient  period  the  The  start-up  can  be  realised  “self-consistently”  (i.e. 

system  is  simulated  in  the  same  way  as  the  first  strategy.  utilising  only  the  energy  produced  in  the  power  module) 
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Fig.  8.  Analysis  of  strategy  no.  3:  stack  temperature  and  maximum  power  furnished  during  the  analysed  transient  period. 
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Fig.  9.  Temperature  and  relative  humidity  start-up  limits  for  strategies  no.  3  and  4.  Above  the  dotted  line,  start-up  can  be  realised  maintaining  humidity 
>50%.  Below  the  solid  line,  start-up  cannot  be  reached  without  humidification  or  without  heating  the  stack. 


in  a  wide  range  of  conditions.  In  spite  of  this  result,  the 
analysis  of  the  trend  of  temperature  and  air  humidity  in  the 
stack  after  the  early  running  stages  shown  in  Fig.  10  indi¬ 
cates  that  the  air  humidifier  is  essential  after  the  start-up. 

The  presence  of  the  air  humidifier  is  also  necessary  to 
avoid  the  membrane  drying  out.  In  fact,  when  drying  occurs 
the  polymeric  film  can  be  broken  and  the  cell  damaged 
irremediably. 

This  work,  however,  shows  the  existence  of  climatic 
conditions  favourable  to  start-up  without  humidifier.  In  fact, 


for  T  =  293  K  and  a  >  50%  (a  typical  weather  condition  of 
southern  Europe),  the  air  humidification  is  not  necessary  for 
the  start-up. 

Contrarily,  for  climatic  conditions  that  can  be  frequent  in 
northern  Europe  (T  =  273  K,  a  <  40%)  the  system  needs 
the  humidifier  and  also  an  air  heater. 

This  consideration  demonstrates  that  the  design  of  the 
power  module  utilising  the  proposed  novel  approach  can  be 
successfully  realised  considering  climatic  conditions  where 
the  application  will  work. 


Fig.  10.  Start-up  operative  limits.  Solid  line  represents  the  external  limit  conditions  for  start-up  in  the  absence  of  air  humidification.  Temperature  and 
humidity  in  the  stack  start  from  the  values  represented  by  the  solid  line  and  follow  the  paths  identified  by  the  arrows.  When  the  arrows  reach  the  solid  line,  the 
power  module  is  not  able  to  sustain  itself:  above  the  broken  line,  obtained  by  matching  the  final  values  of  temperature  with  the  initial  humidity  ones,  the 
system  cannot  operate  without  air  humidification  or  stack  cooling. 
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5.  Conclusions 


Appendix  A 


The  planning  of  a  non-conventional  hybrid  bus  is  an 
engineering  theme  of  great  importance.  The  design  approach 
of  the  European  team  has  shown  that  the  possibility  of 
successfully  constructing  a  working  vehicle  is  dependent 
on  finding  solutions  not  commonly  espoused  in  the  auto¬ 
motive  world,  even  for  accessories. 

In  this  work,  we  have  presented  an  innovative  solution  of  a 
“power  module”  fuel  cell  pack  independent  of  the  vehicle. 
In  particular,  the  start-up  of  this  modular  solution  is  strongly 
conditioned  by  the  air  compression  system,  which  has  been 
shown  to  be  critical  when  its  motor  is  fed  directly  by  the 
electrochemical  generator. 

In  addition,  the  choice  of  this  device  is  not  simple,  given 
the  restrictions  on  the  particular  characteristics  of  air  purity 
the  stacks  require. 

The  course  followed,  however,  permitted  us  to  define  a 
procedure  that  rationalised  the  analysis  of  the  power  module 
under  study,  with  particular  attention  to  the  start-up  pro¬ 
blems  of  the  system  and  the  need  to  implement  some 
auxiliaries  (i.e.  air  humidifier  and  air  heater). 
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Following  Springer  et  al.  [14],  membrane  overpotential  is 
hardly  conditioned  by  the  membrane  hydration,  which 
depends  on  the  electro-osmotic  drag  due  to  the  proton 
crossing  the  membrane  in  solvated  form,  water  diffusion 
due  to  the  non-uniform  value  of  X  and  pressure  gradient. 
Eq.  (14)  includes  these  aspects. 

Boundary  conditions  of  Eq.  (14)  consider  the  evaluation 
of  the  surface  humidity  as  a  function  of  the  external  one: 


Z  =  0, 

O 

II 

(A.l) 

where  X  is 

that  expressed  by: 

0 

X  =  y/l  +  \J/2a  +  1 J/3a"  +  \ J/4a 

(A. 2) 

and 

Pm  0 
a  =  — — 

(A. 3) 

T 


The  integration  of  Eq.  (14)  is  difficult  and  this  is  due  to  the 
fact  that  diffusion  has  been  expressed  by  Springer  et  al.  [14] 
as  a  function  of  X\ 


Dx 


2  3 

=  (cpy  +  (p2X  +  T  cp^X  ) 


exp 


<P  5 


X  >  4 


(A.4) 


The  analysis  of  experimental  data  [14]  reported  in  Fig.  11 
justify  considering  DA  constant  for  X  >  4. 

Adopting  this  consideration,  Eq.  (14)  is  analytically  inte¬ 
grated  by  defining: 


a  =  (A. 5) 

neF 

P'  =  Ca  (A. 6) 


Fig.  11.  Measured  intradiffusion  coefficient  (Df)  adapted  from  [16]. 
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(A. 7) 


3nOPc  —  Pa 

S  =  ky  — - 

A'it  1  m 


y 


(A. 8) 
(A. 9) 


Considering  Eqs.  (A.5)-(A.9),  Eq.  (14)  is  written  as: 


dk 

dz 


a- 


1 


(A. 10) 


The  solution  of  Eq.  (A.  10)  is: 

X  =  ^  [1  -  exp  (&)\  +  20  exp  (&)  (A.ll) 

In  addition,  introducing  the  knowledge  of  the  value  of  X  in 
correspondence  to  the  electrode  (i.e.  X  =  X\  for  z  =  m)  as  a 
function  of  the  water  transport  through  the  membrane,  we 
obtain: 


.  _  -  20exp  (£m)]£y 

a[  1  —  exp  (£m)] 

The  value  of  X  is: 


(A. 12) 


X\  -  X0  exp  (£m) 

2  =  — -  1  ~  exp  (&)  +  20  exp  (fr)  (A.13) 

1  -  exp  (cm) 

So,  introducing  Eq.  (A.13)  into  Eq.  (14)  and  integrating 
Eq.  (13)  the  expression  of  the  membrane  Ohmic  resistance  in 
explicit  form  is: 


exp  [y(T)][/‘(^m)]{^m-ln[(y12i-y2)(exp(^m)-l)] 

r.  =  +  In  [f(£,m)g(X q)]} _ 

^[exp  (£m)g(X 0)  -  2^!  +  72] 

(A.  14) 

where: 

rj-i 

l(T)  =  y  (A.15) 

74  T 

/(£,  m)  =  exp  (£m)  -  1  (A.  16) 

g(20)  =  207i  -  72  (A. 17) 


A  brief  parametric  analysis  of  Eq.  (A.  14)  reveals  that 
diffusivity  does  not  strongly  condition  the  Ohmic  resistance 
evaluation.  In  fact,  by  differentiating  Eq.  (A.  14),  for 
T  =  343  K  and  P  =  1.5  x  105  Pa,  we  obtain: 


7^- <  8000  (A.  18) 

dDx 

and  consequently 

Q<no  +  MOO(Dx-Dj  (A. 19) 

where  Q()  represents  the  minimum  value  of  the  resistance 
calculated  by  Eq.  (A.  14). 

In  addition,  the  maximisation  of  Eq.  (A.  19)  with  the  data 
of  Fig.  11,  reveals  that  membrane  Ohmic  resistance  is  not 


very  sensitive  to  the  diffusivity: 

Q  <  Q0  +  8000(0.5  x  1CT10)  «  Q0  +  0.5  x  10“6  (A.20) 

The  difference  between  the  results  based  on  this  simplified 
approach  (i.e.  Eq.  (A.  14))  and  the  one  obtained  by  the 
numerical  integration  of  Eq.  (A.  14)  (i.e.  where  Dx(X ))  is 
<5%. 
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